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Abstract 

This papers deals with overall phase transformation kinetics. The Fokker-Planck type equation 
is derived from the generalized nucleation theory proposed by Binder and Stauffer. Existence of 
the steady state solution is shown by a method based on the mean value theorem of differential 
calculus. From the analysis of asymptotic behavior of the Fokker-Planck type equation it is known 
that the number of clusters having the critical size increases with time in the case of constant 
driving force. On the basis of the present study on overall phase transformation kinetics a simple 
method for analyzing experimental phase transformation curves was proposed. 
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I. INTRODUCTION 



Overall phase transformation kinetics can easily be detected. A number of techniques 
allow in situ detection. Usually the results are plotted as changes in physical property, such 
as the dilatation and the electrical resistance, versus time. Some of those observation can 
be transformed into the volume fraction of newly precipitated phase versus time. 

The progress of isothermal phase transformation kinetics can be conveniently representing 
by plotting the fraction transformed as a function of time, phase transformation curve. The 
Avrami type of equation [|I| has been utilized for the analysis of given phase transformation 
curve. The fraction transformed f{t) is described by a equation with the form 

/(t) = l-exp(-Cr) (1) 

where C and n are constants which are dependent on temperature and a composition of 
material. Eq.(0) is reduced to 
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In 



In 



C + nlnt (2) 



We assume linear relationship between \nt and ln[ln(l/(l — f(t)))] and then determine a 
numerical exponent n by the slope of the line. The value of n varies from 1 to 4 and it 
depends on dominating mechanism of phase transformation. However applicability of the 
Avrami equation is very limited. 

The overall phase transformation kinetics is known to be descibed by the Kolmogorov- 
Johnson-Mehl- Avrami (KJMA) equation |[^^[Q(see Appnedix A). 

ft 



fit) = 1 - exp 



J{t')V{t,t')dt' 



(3) 



where f{t) is the fraction transformed, J{t) is the nucleation rate at time t and V{t,t') is 
the volume at t of a stable nucleus formed at time t'. The Avrami equation is the special 
case of the KJMA equation in which we assume 

J{t) = Co5{t) 

v{t,t') = Cie 

where S{x) is the Dirac delta function and Cq and Ci are constants. In the above case, we 
assume all nuclei are formed at time zero and the stable nuclei grow with a constant speed. 



In this paper we will deal with a simple model for overall phase transformation kinet- 
ics. First the Fokker-Planck type equation will be derived from the generalized nucleation 
theory proposed by Binder and Stauffer^. Asymptotic behavior of the Fokker-Planck type 
equation will be investigated by using the general mean value theorem of differential calcu- 
lus. On the basis of the KJMA equation a simple method for analyzing experimental phase 
transformation curve will be described. 



II. ASYMPTOTIC BEHAVIOR OF THE FOKKER-PLANCK TYPE EQUATION 

A. Cluster dynamics approach toward generalized nucleation theory 

Binder and Stauffer]^, ^ proposed a dynamical theory of cluster formation in which 
dissociation-coagulation mechanisms are caused by " cluster reactions" . It deals with a wide 
range of cluster sizes, which include both small and large clusters. This theory is considered 
to be a natural extension of the Becker-Doring theory 0]. Binder and Stauffer considered 
the time evolution of the cluster size distribution n{l,t) in terms of a cluster coagulation or 
a cluster splitting mechanism. If we ignore cluster diffusion, we can write an equation for 
evolution of the number of clusters of size 1, n{l,t), 

i'=i i'=i 
1 

+ -J2Ci-,i'n{l',t)n{l-l',t) 
i'=i 

oo 

- Y,CiMl,t)n{l',t) (4) 
i'=i 

The first term on the right hand size of Eq.(^ accounts for an increase in the number of 
/-clusters due to the splitting reaction [/ -|- /'] [I] + [/']. This is assumed to be proportional 
to the number of clusters of size (/ + /') which are present, where Si^in' is the rate constant. 
The second term represents a decrease in the number of /-clusters caused by the splitting 
reaction [/] -—>■ [I — I'] + [/']. The factor 1/2 in the second and the third terms accounts for 
overcounting pairs in the summation. The third term represents an increase in the number 
of / clusters by coagulation of clusters [/ — /'] + [/'] [I]. This is a reverse process of that 
represented by the second term where Ci^ini is the rate constant. The fourth term represent 
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a decrease of /-clusters by coagulation [/] + [/'] [I + 1]. The next step is to assume the 
detailed balance condition. The rates 5* and C can be replaced by a single reaction rate as 



Si+, + 1') = Ci iMi)Mi') = w{i, I') 

where no{l) is the equilibrium cluster distribution function and given by 

e{l) 



(5) 



no{l) = rii exp 



(6) 



with 



(7) 



A function e{l) in Eq.(P) is the free energy of formation of a cluster of size I which contains 
a bulk part of magnitude 5/1 and a surface term with a coefficient 7 that is proportional to 
the surface energy. Substitution of Eq. (KI) into Eq. (0) yields 



dn{l, t) 
dt 



v=i 
i-i 



no{l + V) 



Ml) 

+ ^ n,{l')n,{l-l') 

n{l,t)n{l',t) 



i'=i 



i'=i 



no{l)no{l') 



(8) 



Let us consider the case I » I'. If we assume that dissociation-coagulation mechanisms are 
controlled by small clusters such that /' < [1/2], Equation (||) is rewritten as 



dn{l, t) 
di 



[1/2] 

i'=i 

[1/2] 



no{l + n 
noil) 



[1/2] 



+ ^w^(i-r.n "";;'"';r';f 

^ no{l')no{l - I') 

n{l,t)n{l\t) 



i'=i 

[1/2] 



i'=i 



no(/)no(/') 



(9) 
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We assume that the concentration of small clusters are that of the supersaturated one-phase 
equilibrium. The deviation from equilibrium of t)/no(/') is neglected in the nonlinear 
terms, i.e. n{l',t)/no{l') = 1. Equation (y) is simplified as: 



dn{l, t) 
dt 



[1/2] r 

E 

l'=l 



wiu') 



+ [W{l-l',l') -w{l,l')] 



n{l + l',t) n{l,t) 
no{l + I') ~ no{l) 

n{l-l',t) n{l,t) 



_ noil - I') no{l) \_ 



By expanding the terms n{l + t)/ 71^(1 + /') and W{1 — I', I') about /, we obtain 



n{l±r,t) _ n{l,t) d 
no{l±l',t) ~ no{l) m 



+ 



n{l,t) 



W{l-l',l') = W{lJ')-l 



,dw{i,n 



dl 



Substitution of Eqs. (pT]) and ([1^) into Eq.(|TOD yields a Fokker-Planck type equation 



dn{l,t) 
dt 



d_ 
dl 



d 



n{l,t) 



(10) 



(11) 



(12) 



(13) 



where we have defined a cluster reaction rate a{l) in terms of W{1, V) describing reactions 



[1/2] 



noil) f^. 



(14) 



Thus the Fokker-Planck type equation can be derived in the general case that we allow 
growth and shrinking not only by evaporation/condensation of monomers but also by small 
clusters of size This indicates that we have formally the same expressions of steady state 
cluster concentration and time-dependent nucleation rate as those given by the Becker- 
Doring theory. 

If we consider the contributions of larger clusters, i.e. /' > [1/2], which is neglected the 
above treatment, Eq.(|l^) is rewritten as 



dn{l,t) 
dt 



d_ 
di 



a{l)no{l) 



d_ 

di 



n{l,t) 
no{l) 



+ 



W{l-l',l 



, , n{l, t)n{l — I', t) , 



%(/)%(/ - /') 



dV 



, ^ ' ^ no(/)no(/') 



(15) 



The first term in Eq.(|T5|) is the Fokker-Planck type equation. The last two terms cor- 
responds to standard coagulation equation called Smoluchowski equation. Thus Eq.(|l5]) is 
known to be a combination of two classic equations. 



B. Steady state solution of the Fokker-Planck type equation 

First, we will show existence of a steady state solution of the Fokker-Planck type equation. 
Equation (|13D is rewritten by introducing a new variable 

n{l,t) - ns{l) 



u{l,t) 



as 



no{l) 



du d 



noil) 



(16) 



du 
'dl 



(17) 



where ns{l) is a time independent distribution function. The boundary conditions for Eq(|T^ 



are 



u{l, t) = 0, 1^0, u{l, t) = 0, / ^ oo 



(18) 



Unless u{l,t) = 0, a solution u{l,t) of Eq. (p^) satisfies the above boundary condition has 
at least one positive maximum value and one negative minimal value (see Fig.|I|). 

Consider temporal behavior of a solution along peak top position of u which satisfies 
du/dl = 0. For / = h 



du 



dp 



A function F{t,l,du/dl,d'^u/dl'^,u) is defined as 



du d^u 



' ' dV dP' 

This function satisfies the following condition: 



u 



d_ 
dl 



< 



(19) 



(20) 



F(t,/i,0,0,M) = 



(21) 



We utilize a fundamental relation between the difference quotient and derivative known 
as the mean value theorem of differential calculus p|. Let f{x) is a continuous function 
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nil 




FIG. 1: Schematic illustration of trajectories of a peak top and a bottom of a valley of a function 
u{l,t) 

on the closed interval [a,b\. We assume that the derivative exists everywhere in the closed 
interval[a, h]. There exits an intermediate value ^ in the interval such that 

/(&) - /(«) 



b — a 



no 



(22) 



This statement is called the mean value theorem of differential calculus. Prom this ordinary 
mean value theory we can derive the general mean value theorem of differential calculus. 

Wc assume u{x) is monotonic in the interval [x,y]. The compound function f{u{x)) for x is 
defined for u{x) in the interval [■u(x), We find that there exist a function /'(^) for an 

intermediate value ^ in the open interval {u{x),u{y)) such that 



f{u{y))-f{u{x))=r{i)[u{y)-u{x)] 



(23) 



We obtain the following equation for an intermediate value C, in the open interval 
{{d''u/df),Q) 

( du\ ^ f , du d'^u \ , 
dF{t,hAC,u) 

dp J d{d^u/dP) ^ ' 

This equation indicates that the function u along the peak top position decreases with time 
and approaches the constant valueci . By repeating the above procedure for a solution along 







vO, t) 




FIG. 2: Schematic illustration of a trajectory of a bottom of a valley of a function v{l,t) 



the valley bottom position, it is known that the negative minimum value increases with time 
and approaches the constant value C2- 

As the absolute values of u along peak top and valley bottom positions decreases with 
time, the function u is considered to be time independent at longer times. We will show 
that the constant value Ci is equal to C2- By differentiating the both members of Eq.(|l] 
with respect to I we obtain 

du 



a{l)no{l) 



dldt 



aW||ln«„(/)|| 



(25) 



By introducing a new variable v, 



V = a(/)no(/) 



du 
'dl 



Eq.(P5|) is rewritten as 



dv 
di 



(26) 



V is negative at any point within the open distance (Zi, I2) and v{li) = v{l2) = 0(see Fig.||). 
Let us consider the trajectory of v{l,t) along the bottom of the valley, /3(t) where 

dv{i,t) ^ dMi,t) 



dl 



From Eqs. dD, (^), it follows that 



dv{l, t) 
dt 



> 



> 



(27) 



(28) 
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This indicates that the value of v becomes at longer times. Then it is known that the 
value of u approaches and we have the steady state solution. 

Becker and Doring used the following choice of boundary conditions to obtain time- 
independent solution : 



noil) 
noil) 

The steady state solution of Eq. (|T3|) is 



0, if l^oo. (29) 
^ 1, if 1^0. (30) 



Jo a{l)no{l) 

The range of integration in Eq.(|3T|) include the peak at I = Ic where e(/) is maximum and the 
integral can be evaluated by the saddle point approximation. The steady state nucleation 
rate, tT^q, can be written as follows: 



dhil)/dn=iX^' , ^ r eik 



Js = ailc) ( 2iTkBT ~ ' ^^^^""^ ^ ^oexp 



(32) 



where Jq is the so called nucleation prefactor. The quantity e(/c) is an activation energy 
(the energy for forming critical nucleus). The factor i—d'^eil)/dP\i=ij2nkBT)^^'^ in the right 
hand side of Eq.(p2D is known as the Zeldvich Factor, Z. 

The above mentioned steady state nucleation rate does not provide any information on 
the momentary cluster distribution nor on the nucleation rate prior to reaching the steady 
state condition, i.e time dependent nucleation rate. It was found |TD| that the nucleation 
rate. Jit), approaches the steady state 

J(t) = J,[l -exp(-At)] (33) 

where 6.3a(/c)2'^ < A < 12.0a(/c)^^- This result agrees with numerical calculation by 
Kanne-Dannetschek and Stauffer where A = 6.4a(/c)2'^ 



C. Nucleation kinetics in materials with a time-independent driving force 

As shown previous section, the nucleation rate approaches a constant value, steady state 
nucleation rate, because the driving force of nucleation in an infinite system is time de- 
pendent. Here let us consider the case that the driving force is constant during phase 
transformation. 
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The Fokker-Planck type equation([T5|) is rewritten as 



dn{l,t) dMl^ 

= — — \- 



dt 



+ 



dP 
1 r dajl) 



+ a{l) 



da{l) ail) deil) 
dl ^ keT dl 



dn{l,t) 
dl 



n{l,t) 



(34) 



keT \ dl ' dP 

If the driving force of nucleation is constant the value of e{l) becomes maximum at Z = Z 
then we have 



dl ~ ' dl 



0, 



(35) 



From the assumption that the driving force is time-independent, it follows that the minimum 
value of n{l, t) is given at / = Ic- 



dl 

From Eq.(0) it follows furthermore that 

d^ejl) 
dP 

Because very large clusters do not exist 



dn{lc,t) ^ ^ d'^njlct) ^ ^ 



dp 



0, I ^ oo 



(36) 



(37) 



n{l, t) ^ 0, I ^ oo 



(38) 



we obtain 



dMl,t) 
dp 



0, / ^ oo 



(39) 



Now we consider behavior of n{l, t) along the trajectory peak top position of n{l, t), I = lc{t). 
A fuction G{tJ,n,dn/dl,d'^n/dP) is difined as 



. dn d'^n 

''''~dl'~dP 



d_ 

dl 



a{l)no{l) 



d_ 

dl 



n{l, t 
no{l) 



(40) 



The fuction G satisfies 



G(t,/e,n(/e,t),0,0) = 



(41) 



Applying the mean value theorem of differential calculus for compound functions to the 
function G, we obtain the following equation for an intermediate value ( in the open interval 
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FIG. 3: Schematic illustration of a trajectory of a bottom of a valley of the cluster size distribution 
fuction n{l,t) in the case of constant driving force 



From Eg. (^2]) it follows that the number of stable clusters increases with time in the case 
that the driving force is time- independent. (see FigJ^). 

III. OVERALL PHASE TRANSFORMATION KINETICS 

Here we propose a method for analyzing given phase transformation curve with use of 
the KJMA equation. This method is based on the following assumptions. 
(l)The time dependent nucleation rate, J if), is described by 



(0,^2 



n/df) 




G(t,/e,n(/„t),0,0) 
dP d{dydP) 



> 



(42) 



J(t) = JJl - exp(-At)] 



where Jg and A are constants. 
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(2) The volume at t of stable nuclei formed at time t', V{t^ t') is a function of t — t', that is 
V{t,t') = V{t-t'). 

As shown in the previous section, the assumption (1) is reasonable one and the assumption 
(2) is acceptable unless the supersaturation of metastable solid solution changes rapidly. 

First we determine the variation of the extended volume, V^x = Jq ■J{t')V{t,t')dt' , with 
time by given phase transformation curve. 

1 



In 



[I- fit) 



From the assumptions (1),(2), it follows that 



Vexit) = J, 



[ V{t-t')[l-exp{-At')]dt' 
Jo 



(43) 



(44) 



By perfoming change of variable as t" = t — t' we obtain 



J. 



V{t")dt" 



- exp(-At) f V{t")exp{At")dt" 
Jo 



(45) 



By differentiating both members of Eg. (^51) with respect to t and with use of the relation 
V{Q) = we have 

exp(y4t) dVex 



AJs dt 

Further differentiation of both members of Eq 



V{t") exp{At")df 



with respect to t yields 



V{t) 



J.A 



A- 



dKr. 



+ 



d'^V.s^ 



(46) 



(47) 



dt ' dt'^ 

where the constant Js is an adjustable parameter which can be determined by given phase 
transformation curve and the parameter A can be estimated by the numerical calculation 
of the Fokker-Planck type equation^ to be 6.4a(/c)Z^. Thus we can determine the change 
in the volume of stable nucleus with time by given phase transformation curve as shown in 
Fig.|. 

Next we consider a method to estimate the size distribution function by given phase 
transformation curve. Let A/(to) to be a volume of nuclei formed at a interval of time 
between to and to + ^t. The contribution of nuclei formed after a time t" to the fraction 
transformed at time t, /'(t,t"), is given by 



1 — exp 



J{t')V{t,t')dt' 



(48) 
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FIG. 4: Estimation of the volume of stable nucleus V{t) by a phase transformation curve, f{t) 
versus time 

Then we obtain A/ (to) as 

A/ (to) = f'{t,to)-f'{t,to + At) 
— exp 



As At ^ 



— exp 



-Js f V{t - t')[l - ex.p{-At')]dt' 

Jto+At 

-Js [ V{t-t')[l-e^p{-At')]dt' 

Jto 



A/(to) d 



lim 
At^o At 



dto 



exp 



Js I V{t-t') 

Jto 



[1 - exp(-At')]cit' 



(49) 



(50) 
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FIG. 5: The size distribution function estimated by the phase transformation curve shown in Fig. 



Equation (EDf) yields 



At~*0 At 



J,V{t-to)[l~exp{-Ato)] 



exp 



Js [ V{t - t')[l - exp{-At')]dt' 

J to 



(51) 



At a interval of time between to and to + At, Js[l — exp{—Ato)] nuclei formed per unit 
volume. As the fraction transformed at the time to is /(to), the number of nuclei is given 
by Js[l — f{to)][l — exp{—Ato)]. Hence the average volume at time t of nuclei formed at a 
interval of time between to and to + is given by 



Vit - to) 



exp 



-Js V{t-t')[l-exp{-At')]dt' 

J to 



(52) 



1 - /(to) 

The fuction V{t — t') in Eq.(|S2) can be determined by given phase transformation curve. 



With use of Eq.(^) we can estimate the size distribution function of nuclei. Figure ^ shows 
the size distribution of nuclei estimated by the phase transformation curve in Fig. 
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IV. CONCLUSION 



In this paper we attempted to analyze asymptotic behavior of the Fokker-Planck type 
equation which is derived from the generahzed nucleation theory by Binder and Stauffer. 
From the analysis it follows that the number of cluster having the critical size increases with 
time if the driving force is constant during phase transformation. However, the chemical 
driving force which mainly controls diffusional phase transformation decreases with time. 
Therefore, the number of stable clusters increases with time in the case that an external 
driving force is applied to compensate the decrease in the chemical driving force. One of the 
most effective way is to give the large amount of plastic deformation to metastable phase 
at higher supersaturation. This may be helpful to develop a new method for producing 
structural materials. 

Next on the basis of a model for overall phase transformation kinetics, we propose a 
simple method for analyzing experimental phase transformation curves. This method is 
based on the Kolmogorov-Johnson-Mehl-Avrami equation. Time-dependent nucleation rate 
is determined by the analysis of asymptotic behavior of nucleation. A function for describing 
the diffusional growth can be estimated by the phase transformation curve. Furthermore 
the size distribution function of newly precipitated phase can be predicted with use of the 
present analysis. 

APPENDIX A: GENERALIZED KOLMOGOROV-JOHNSON-MEHL-AVRAMI 
EQUATION 

We derive a formula for the time evolution of the volume fraction of the transformed 
phase at t, f{t)^. The volume fraction of the metastable phase is given by 

git) = 1 - fit) (Ai) 

Take a point A at random. We calculate the rate of the fraction transformed at the point 
A within the period t to t + dr 

dfir) = -dg{r) = -[g{T + dr) - gir)] (A2) 

Let v{P, t) be the growth rate at time t of a droplet which nucleated at a point P at a time 
t' . The distance between the point P and a point Q, which is the intersection of the line AP 
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FIG. 6: Condition for a point A reached by a droplet at moment t' within the period between r 
and T + dr 

and the surface of the droplet (see Fig.|]), at time t" is given by 



PQ = Z{P,t",t') = v{P,t)dt (A3) 



t' 



The surface of a droplet formed at a point P' at moment t' reaches the point concerned 
within the period t to t + dr only if the distance between A and P' satisfies the following 
condition 

I v{P',t)dt <P^< / v{P',t)dt (A4) 
Jt' Jt' 
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Thus the point A is reached by droplets formed at time t' within the period between r and 
T + dr only if nucleation occurs points within a closed domain of volume 



= j J Z{P',T,t')v{P',T)dxdydT 



(A5) 



where D is the surface on which each point is situated at a distance Z{P'{x, y, z),t, t') from 
the point A. The probability of nucleation in this domain within the period t' to t' + dt' is 
equal to 



P^dt' = j j Z{P', T, t')v{P\ T)J{t')dxdydt' 



(A6) 



where J{t') is nucleation rate per unit volume at time t' . Then the probability that droplets 
formed at a randomly selected point during the period of time between t' = and t' = t 
reach the point A is 



Pa= f [ [ Z{P',T,t')v{P',T)J{t')dxdydt' 
Jo Jd j 



(A7) 



The point A transforms to a stable phase from the metastable phase only if the point has 
not transformed previously. Thus, the rate of transformed is 

dfir) = -dg{r) = g{r)PA (A8) 

The solution of Eq.(|A5|) is given by 

\n^ = - [ dr [ [ [ Z{P',T,t')v{P',t)J{t')dxdydt' (A9) 
5'(0) Jo Jo Jd j 

Since the integral of the right hand side of Eq.( |A9|) has a finite value, we can change the 

sequence of integration. Then 

In 44 = -/" J(t')dt' [ [ [ Z(P',T,t')v(P',t)dxdydT (AlO) 
^(0) Jo Jt' JdJ 

The volume at a time t of a droplet which formed at a time t' is 

v{t,t') = J^^J J Z{P',T,t')v{P',t)dxdydT (All) 
At t = 0, g{0) = 1 , we obtain the well known Kolmogorov-Johnson-Mehl-Avrami type 



formula 



fit) = 1 - exp 



- [ J{t')V{t,t')dt 
Jo 



(A12) 
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